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ABSTRACT 

We present high resolution Na I D spectroscopy of the /3Pic disk, and the resonantly scattered sodium 
emission can be traced from less than 30 AU to at least 140 AU from the central star. This atomic gas 
is co-existent with the dust particles, suggestive of a common origin or source. The disk rotates towards 
us in the south-west and away from us in the north-east. The velocity pattern of the gas finally provides 
direct evidence that the faint linear feature seen in images of the star is a circumstellar disk in Keplerian 
rotation. From modelling the spatial distribution of the Na I line profiles we determine the effective 
dynamical mass to be (1.40 ± 0.05) M©, which is smaller than the stellar mass, 1.75 M©. We ascribe this 
difference to the gravity opposing radiation pressure in the Na I lines. We argue that this is consistent 
with the fact that Na is nearly completely ionised throughout the disk (Na I/Na < 10~ 4 ). The total 
column density of sodium gas is A(Na)=10 15 cm~ 2 . 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: individual 
08 Pic) 



1. INTRODUCTION 

The /3Pic disk has been the subject of intense studies 
ever since it was shown to have excess emission in the ther- 
mal infrared by the IRAS mission (Aumann 1985). The 
dust particles giving rise to the infrared emission also scat- 
ter the stellar light and coronographic techniques revealed 
an elongated structure (Smith and Terrile 1984) that re- 
cently has been extremely well characterised by means of 
HST observations (Heap et al. 2000). Although the gen- 
eral impression of these images is a high degree of smooth- 
ness and symmetry, Heap et al. (2000) show that the peak 
of the elongated scattered emission has wiggles that sug- 
gests that the inner part of the disk is tilted relative to 
the outer part. There is also evidence for asymmetries 
and structures in the radial light distribution of this inner 
disk. Asymmetries have also been observed in the mid-IR 
(Pantin, Lagage and Artymowicz 1997). 

These asymmetries and the remarkable discovery 
(Vidal-Madjar et al. 1994; Lagrange, Backman and Arty- 
mowicz 2000) of comet-like bodies occasionally occulting 
the star and giving rise to high velocity absorption lines 
reveal dynamical activity that generally is interpreted as 
evidence for planetary disturbances (Artymowicz 1997; 
Augereau et al. 2001). 

When it comes to the gas component in the disk, the 
situation is at present confusing. Absorption spectroscopy 
in the optical and UV shows a stable component at the 
same radial velocity as the star (Vidal-Madjar et al. 1986; 
Lagrange et al. 1998). One would expect radial expansion 
of this gas, caused by the radiation pressure of the star 
and as a possible explanation a ring of atomic hydrogen 
has been proposed (Lagrange et al. 1998). If present, such 
a ring must be close to the star as widespread H I is not 
detected at 21 cm (Freudling et al. 1995). CO has been de- 
tected by means of UV spectroscopy (Roberge, Feldman 
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and Lagrange 2000) but sensitive searches for CO emis- 
sion in the radio region have failed to detect this molecule 
(Liseau and Artymowicz 1998). Molecular hydrogen is not 
detected in UV absorption spectra (Lecavelier des Etangs 
et al. 2001) but on the other hand, ISO/SWS spectroscopy 
(Thi et al. 2001) indicates H2 emission lines that, if real 
(the S/N is low), require large amounts of molecular hy- 
drogen in a spatial distribution that has a void in the line 
of sight towards the star. 

So far it has not been possible to directly probe the ve- 
locity pattern along the disk, lacking the detection of spa- 
tially resolved gas emission (except a possible detection of 
Fe II emission very close to the star, Lecavelier des Etangs, 
Hobbs and Vidal-Madjar (2000)). In the present paper we 
report the detection of the resonance Na I doublet lines at 
5990 A and 5996 A. 

2. OBSERVATIONS AND RESULTS 

/3Pic was observed on January 2-4, 2001, with the 
Echelle spectrograph EMMI on the 3.5 m NTT (New Tech- 
nology Telescope) at ESO-La Silla, Chile. For the Na I 
D1/D2 observations, we used a 5arcmin long slit and 
therefore omitted the cross-disperser. Instead, we inserted 
an order-sorting interference filter, which was centred on 
5890 A and 60 A wide. We have measured the filter pro- 
file and are in control of the leakage from the neighbouring 
orders (less than 3% and 4% respectively). With the wave- 
length resolution of 6 • 10 4 (At> = 5 km s _1 ) this instrumen- 
tal set-up can be viewed as some sort of a 'spectroscopic 
coronograph', where the high dispersion is used to prevent 
the saturation of the detector (CCD) by the very bright 
stellar point source. Still, the overheads of the observing 
time were dominated by detector readout times, since in- 
tegration times of individual observations were typically 
only 300 s. The total integration time was chosen such 
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that the co-added spectra should permit the detection of 
the faint light scattered off the disk by the dust, which 
was also achieved. Parallel with the disk (at position an- 
gle 30°75, Kalas and Jewitt (1995)), this time was 2.5 hr, 
whereas we spent in total 36 min for observations with the 
slit oriented perpendicular to the disk. The seeing values 
determined from these observations are 1"2, whereas the 
width of the slit was kept at l'.'O in order to sample (on the 
CCD) the spectral resolution at the Nyquist frequency. 

In addition to the two telluric emission lines, the paral- 
lel observations also show spatially extended Na I D emis- 
sion, both lines consistently redshifted to the heliocentric 
velocity of /3Pic (+21 kms^ 1 ). These observations are dis- 
played in Fig. 1. In Fig. 2, the distribution of the Na I 
emission along the slit, on either side of the stellar spec- 
trum, is compared to the corresponding distribution of the 
dust. Evidently, both gas and dust show similar light dis- 
tributions, suggestive of the co-existence of these species 
in the (3 Pic disk. 

3. DISCUSSION 

3.1. The rotation of the disk 

We first address the question of Keplerian rotation. Jus- 
tified by the indications that there is no significant outward 
motion, we simply assume circular rotation. The distance 
to the star is 19.3 pc (Crifo et al. 1997). In order to in- 
terpret the velocity pattern, we must first find a radial 
density distribution of the Na gas component that results 
in the observed light distribution. For the moment, we 
ignore the difference between the two sides (Fig. 2) and 
use the average light distribution. Adopting a functional 
relation that is suitable for a broken power law: 

n(Nal) oc [(r/a) 2b + (r/a) 2c ] ~ h 

where n(Na I) is the volume density of Na I at the radial 
distance r (within the 1" slit), we get a quite satisfactory 
fit to the observed projected light distribution with the 
following numerical values: a — 89.5 AU, b — —1.28 and 
c = 3.44. The distribution function is shown in Fig. 3 and 
the fit to the observations in Fig. 4. 

Our observations are consistent with circular Keplerian 
rotation of the gas (see Fig. 5), were we have taken the 
measured seeing and spectral instrumental profile (using 
telluric absorption lines) into account. The best-fit mass, 
in a x 2 -sense, equals (1.40 ± 0.05) M©. This model mass 
is less than that deduced from the stellar position in the 
HR-diagram, viz. 1.75 M© (e.g. Crifo et al. (1997)). It 
is possible that this difference is due to radiation pressure 
(see Sect. 3.3). 

3.2. The light budget 

The origin of the "stable gas" absorption has been ex- 
tensively discussed and there are arguments for a location 
very close to the star (Lagrange et al. 1998). If, on the 
other hand, the sodium emission far out in the disk is 
caused by resonance scattering of stellar radiation then, 
at least partly, the "stable gas" absorption must originate 
throughout the disk along the line of sight. We measure 
an Na ID2 equivalent width in absorption of 9.4 mA, im- 
plying a column density of iV(Nal) = 7 • 10 10 cm~ 2 , in 
agreement with previous observations by Vidal-Madjar et 



al. (1986). In emission, a total equivalent width of 0.72 mA 
is obtained. This latter number is of course a lower limit 
as it does not include any line emission originating out- 
side the one arcsecond wide slit of the spectrograph. It 
means that the disk, as seen from the star, must occupy 
at least a latitudinal angle of 8°8. We start to detect the 
line emission at a distance of 30 AU and at this distance 
the required thickness of the disk would be at least 4.6 AU. 
As the slit width covers 19 AU these estimates do not lead 
to any contradiction regarding the light budget under the 
assumption of resonance scattering. On the other hand, 
there is not much margin for the proposed dense H I ring 
(Lagrange et al. 1998) to contribute to the stable Na I 
absorption lines. 

3.3. The radiation pressure 

How, then, can we understand why the radiation pres- 
sure does not quickly accelerate the Na I outwards in the 
same sense as it does, for instance, in comets (cf. Cre- 
moncse et al. (1997))? We first note that the momentum 
transfer caused by the resonance scattering vastly exceeds 
the gravity. Assuming a stellar mass of 1.75 M Q we find 
that the force caused by the radiation pressure exceeds the 
gravitational force by a factor of 300. A possible explana- 
tion why we do not observe a radial velocity component 
is the presence of a relatively dense gas component. How- 
ever, this would require very large amounts of gas far out 
in the disk. As an example we consider the distance of 
100 AU. Our observations exclude a radial velocity com- 
ponent exceeding 5kms _1 and assuming that the main 
gas component is atomic hydrogen we find that the gas 
density must be at least 10 5 cm~ 3 . As the thickness (or 
scale height) of the gas component of the disk is unknown, 
it is hard to judge if this minimum density is in conflict 
with other observations. 

Even though we cannot at present rule out this expla- 
nation for the balancing of the radiation pressure, we pro- 
pose an alternative: Due to the low ionisation potential 
of Na I (5.1 eV), one would expect the stellar UV radia- 
tion to provide a high degree of ionisation. As Na II lacks 
strong transitions within the range of the stellar spectrum, 
the radiation pressure on the ionised sodium gas does not 
counterweigh the gravitation. Thus, if the gas density is 
low, a sodium atom will quickly take up speed outwards. 
But it will stay neutral just for a short time and then, 
as singly ionised, it will have a long time to adapt to the 
motion of the main gas component. 

To quantitatively test this idea we estimate the degree 
of ionisation. As we can exclude higher ionisation stages 
(the ionisation potential is 47 eV for Na II) the equation 
of ionisation equilibrium reads 

rc(Nall) _ w(r)r 
"(Nal) ~ atot [T(r)] n c (r) 
where u> is the dilution factor of the radiation density at 
the distance r from the star compared to that at the stel- 
lar surface, T is the ionisation rate at the stellar surface, 
a t ot is the total recombination coefficient and n e the elec- 
tron density. Using a model atmosphere from Allard et al. 
(2000) and cross sections from Cunto et al. (1993) we es- 
timate Lj = 205 s _1 . Assuming 100 K as a typical electron 
temperature in the disk we get a to t — 3.85 • 10~ 12 enr^s" 1 
(Verner and Ferland 1996). The electron density remains 
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to be estimated. The stellar far-UV continuum does not 
suffice to ionise H, He (if at all present), C, N and O. As- 
suming solar abundances (Holwcger and Rentzsch-Holm 
1995) for elements with ionisation potentials < 8.3 eV 
(longward of 1500A) we estimate n c ~ 50n(NaII). If 
we, finally, assume a constant thickness of the disk (the 
half power width of the dust disk is in fact close to con- 
stant from 30 AU and outwards, see Heap et al. (2000)) 
we can use the radial distribution derived for Na I in com- 
bination with the column density to derive the degree of 
ionisation as a function of the distance to the star. We find 
that sodium is indeed highly ionised with n(NaI)/n(Na) < 
10 -4 throughout the disk. In Fig. 6, we show the radial 
distribution of n(Na I) and n(Na II). These analytical re- 
sults have been confirmed by self-consistent photoionisa- 
tion computations, which produce the observed line fluxes. 
The details will be given in a forthcoming article. 

We conclude that the radiation pressure on sodium, av- 
eraged over time, would be small. We must, however, 
also estimate what speed a neutral atom would typically 
achieve until it becomes ionised. Both the UV contin- 
uum and the sodium D lines are optically thin throughout 
the disk and thus both ionisation rate and the resonance 
scattering rate roughly scale as the inverse square of the 
distance to the star. Therefore, the number of scattered 
photons per period of neutral state is constant and we find 
that typically 1.5 • 10 4 scattering events will occur until an 
atom becomes ionised. The scattered radiation is essen- 
tially isotropic and the transferred momentum would cause 
an outward velocity of merely 0.4 km s _1 . After that, there 
will be a long period of time during which the sodium ion 
will interact with the main gas components and conform 
to the general velocity pattern, be it basically circular Ke- 
plcrian rotation or not. This also means that the radiation 
pressure may cause a slow, stepwise motion of each sodium 
atom outwards resulting in a net flow from the inner to the 
outer parts of the disk. 



3.4. The amount of atomic gas in the disk 

From the radial distribution of Na I we were able to de- 
rive the radial distribution of Na II (Fig. 6). The number 
density is typically n(Na)=lcm~ 3 and the column den- 
sity from 30 to 140 AU is A(Na)=10 15 cm- 2 . If the Na/H 
ratio were solar, this would indicate a hydrogen column 
density of iV(H) = 5 • 10 20 cm -2 . Further discussion of 
these aspects will be deferred to a forthcoming paper. 

4. CONCLUSIONS 

We have observed Na I Dl and D2 emission along the 
disk of /3Pic from a projected distance of 30 to 140 AU. 
The velocity pattern mimics circular Keplerian rotation, 
but the deduced mass is somewhat lower than that ex- 
pected for an A5V star (1.4 compared to 1.75 M Q ). This 
difference is probably due to radiation pressure that to 
some extent counterbalances gravity. We find that the 
sodium gas is mainly ionised, with Na II/Na I around 10 4 
and for this reason, the radiation pressure does not accel- 
erate the sodium gas component to high velocities. This is 
simply because the number of scattering events during a 
neutral period of a sodium atom only would suffice to give 
an outward velocity of 0.4 km s -1 and then, during the 
typically 10 4 times longer period of singly ionised state, 
the radiation pressure is negligible. So, there is plenty of 
time for the sodium gas to conform, through gas friction, 
to the motion of the main gas component (i.e. H, O, C 
and N) that has only a weak direct interaction with the 
radiation field. 

It is, finally, clear that the observations presented in the 
present Letter only mark a starting point of spectroscopic 
investigations as both the spectral and the spatial reso- 
lutions as well as the straylight level can be significantly 
improved upon. In addition, other lines like the Ca II 
H & K doublet and a number of resonance lines in the UV, 
accessible from e.g. the HST, will probably be detected in 
the disk. 
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Fig. 1. — The spectral region around the Na I D lines toward /3Pictoris. The dispersion is along the horizontal axis, with wavelengths 
increasing to the right, and the scale is 0.035A per pixel (1.8kms — ipxl" 1 ). The spatial dimension is along the vertical axis, with the 1" 
wide and 300" long slit oriented along the position angle 30°75. The scale is 0.27arcsec per pixel (5 AU pxl - 1 ), with positive values toward 
the south-west and negative values toward the north-east. The (deliberately very much reduced) stellar continuum at the centre of the figure 
is covered with telluric absorption lines, whereas the terrestrial ionospheric Na I D1/D2 emission lines extend over the entire vertical space. 
The Na I line emission from the /3Pic disk is seen redshifted with respect to the sky lines and with the intensity ratio 1:2. On either side of 
the star, the disk lines are either blue- or redshifted with respect to the systemic rest wavelength, represented by the disk gas absorption seen 
against the stellar continuum. 
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Fig. 2. — The radial distribution, for the two sides of the disk, of the Na I D2 line emission with error bars compared to that of the scattered 
stellar radiation caused by the dust (from Heap ct al. (2000)). The flux scale is arbitrary and adjusted to facilitate the comparison. 
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Fig. 3.— 



The derived radial distribution for Na I (normalised at the peak). 
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Fig. 4. — The projected radial distribution of the Na I D2 emission of the model (smooth line) compared to the observed (average of the 
SW and NE sides). 




Fig. 5. — Comparison of the observed Na I D2 line with the results of theoretical model calculations. To the left, the model is displayed 
at spatial and spectral resolutions 10 times better than the observed. In the middle we show the model degraded to the quality of the 
observations (including noise). As is seen, it compares well with the observations (to the right). 
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Fig. 6. — The derived number density (cm 3 ) for Na I and Na II. Only a small fraction (< 10 4 ) of the sodium is neutral. 



